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After an X-ray binary experiences a transient jet ejection, it undergoes a phase in which its X-ray 
light curve is dominated, for some time, by thermal emission from an accretion disk surrounding 
the black hole. The accretion physics in the thermal-dominant state is understood better than in 
any other, and it is therefore the best state for comparing observations to theoretical models. Here, 
I present simulations that study the way a thermally-emitting disk might be expected to behave 
immediately after a large-scale, steady jet has been removed from the system in the form of a 
sudden ejection. I simulate the ejection's effect on the disk by allowing the strength of turbulence 
(modeled by the a parameter of Shakura and Sunyaev) to increase rapidly in time, and I show 
how this change can lead to an outburst in an otherwise-steady disk. The motivation for treating 
the jet removal in this way is the fact that many models for jets involve large-scale magnetic 
fields that should inhibit the magnetorotational instability believed to drive turbulence; this should 
naturally lead to a rapid increase in turbulence when the magnetic field is ejected from the system 
or otherwise destroyed during the ejection event. I show how the timescale and luminosity of the 
outburst can be controlled by the manner in which a is allowed to change, and I briefly discuss 
ways in which these simulations can be compared to observations of X-ray binaries, in particular 
GRS 1915+105, which shows the most complex and variable behavior of any black hole system 
in outburst. 
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1. Introduction 

Thanks to hard work by many researchers, much is now known observationally about the be- 
havior of X-ray binary systems as they undergo an X-ray outburst coupled with a transient radio jet 
ejection (e.g., [§, 12[ 3, 16, 21, |6||). The evolution during the outburst can be conveniently plotted 



on a hardness-intensity diagram, showing that the ejections occur as the accretion disk moves from 
a bright hard to a bright soft X-ray state and as a steady jet structure becomes unobservable at radio 



wavelengths [12] (although there is a large amount of variation in the detailed tracks that systems 



take on this diagram; see, e.g., [26]). 

The soft state, in which thermal emission from the accretion disk is dominant (hence the 
alternative name "thermal-dominant state") and there is no observational evidence for a jet, is the 
state that is by far the best understood physically. Going back to the original work of Shakura and 



Sunyaev [ )28[ ] (who explained the basic structure of a thin disk), Abramowicz and collaborators [jl|] 
(who extended the theory to include moderately thick disks and showed the possibility of limit- 
cycle oscillations in the inner disk associated with the radiation pressure instability) and many 
others, the basic physics of the disk during this state is reasonably well-known. In trying to compare 
theory to observations, the thermal-dominant state provides the best hope for a stringent test. 

Given the observational picture that has been built up, it is worthwhile to ask the following 
previously-unasked question: What would we expect to happen to a thermal-dominant accretion 
disk immediately after a large-scale, steady jet has been ejected from its environment? 

GRS 1915+105 is the X-ray binary which is best suited to answering this kind of question. It 
has essentially spent all of its time since its discovery (nearly fifteen years ago) in an "active" state, 
cycling right around the region of the hardness-intensity diagram where transient jets form (e.g., 
JT2|]). What makes GRS 1915+105 unique is likely its enormous accretion disk, which is much 
bigger than that of any other similar system and therefore effectively gives the black hole's strong 
gravity a lot of mass that it can continually "play with" [|7, 11, 26]. Multiwavelength observations 
of GRS 1915+105 in the X-rays, infrared and radio show that the system sometimes produces 



repeated ejections on timescales of ~ 30 minutes or faster [|25J, |8|, |23|, [270. An example of this 
behavior is shown in Figure [l| A hard-to-soft state transition occurs around the time of a "spike" 
in the X-ray light curve (at ~ 2, 200 seconds), after which the nonthermal power law component of 
the X-ray spectrum becomes weaker and the thermal disk component becomes stronger. The spike 
also seems to coincide with the beginning of an infrared flare, signaling that an ejection has taken 
place. The timescale over which the state transition happens is extremely rapid, on the order of a 



couple seconds in this case [27] 



Following the ejection, the disk enters a period of repeating oscillations. Numerical simula- 



tions can generally reproduce the basic temporal Q17[ ] and spectral [|33|] properties of these types of 



oscillations using a classic limit-cycle instability. Nonetheless, it is unclear how a limit-cycle insta- 
bility can be used to explain the initial outburst in Figure [j], which is clearly of longer duration and 
a different character than the subsequent oscillations. Nor is it clear what causes the oscillations 
to start and the outburst to occur in the first place. Here, I present a simplified model in which the 
rapid removal of magnetic field associated with a steady jet can actually drive an otherwise-steady 
disk into an outburst whose properties and timescale can be more easily matched to the initial 
post-ejection behavior seen in GRS 1915+105 and other objects. 
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Figure 1: An example of a transient jet ejection in GRS 1915+105 on a timescale of ~ 20 minutes, based 
on observations using the Rossi X-ray Timing Explorer and the Palomar 5 meter telescope [^J]. The top two 
panels show the infrared (dereddened by 3.3 magnitudes) and X-ray light curves, respectively. The bottom 
two panels show the results of X-ray spectral fitting at 4-second resolution. Notice that the infrared flare (sig- 
naling a jet ejection) appears to begin around the time of a "spike" in the X-ray light curve that corresponds 
to a state transition, in which thermal emission from the disk begins to dominate over nonthermal power law 
emission. Notice also that the behavior of the disk immediately after the spike has a different character than 
the subsequent pattern of X-ray oscillations, which show rapid changes in the disk temperature. This figure 



is modified from one originally presented by Rothstein et al. 2005 [ 27 1 



2. Simulations 

I present results from one-dimensional (i.e., axisymmetric and height-integrated) simulations 



of a standard accretion disk around a 14 M e black hole, appropriate for GRS 1915+105 [13, M 



Radiative and advective cooling are included in the equations, turbulence is treated using the dimen- 



sionless a parameter of Shakura and Sunyaev [28], gravity is treated using a pseudo-Newtonian 



potential [gj, and large-scale, ordered magnetic fields are ignored (because they are assumed to be 
removed from the system when the steady jet disappears). More details of the simulations will be 
presented in a forthcoming work, but in general, the physics behind them is similar to that included 
in other one-dimensional disk simulations that have been presented in the literature. 

All of the simulations begin with a gas pressure dominated disk with a dimensionless accretion 
rate M = 0. 1 (measured in units of the Eddington luminosity divided by the speed of light squared). 
For these parameters, the entire disk in the region of the simulation (typically lr g out to a few 
hundred r g , where r g is the gravitational radius) is initially stable. 
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For the purpose of these simulations, the primary effect of the "removal of a jet" is postulated 
to be a rapid increase in a, the parameter that controls the strength of turbulence in the disk. From 
an observational standpoint, this idea is attractive because a is essentially the "clock" in a standard 
disk; the thermal and viscous timescales are inversely proportional to it. Observations such as 
those in Figure [I] clearly show that the disk experiences a state transition on a very rapid timescale, 
and after the transition the typical evolution timescales of the disk are faster. From a theoretical 
standpoint, meanwhile, many models suggest, usually indirectly, that the value of a should be 
different when a jet is and is not present. These models invoke a picture for a steady jet in which 
a large-scale, ordered magnetic field is present that can channel matter and energy away from the 
accretion disk (e.g., [Q, |10[ 31]). If these magnetic fields are stronger than equipartition, the 
magnetorotational instability that is thought to drive accretion is not generally expected to operate 

and the strength of turbulence (parameterized by a) will be weaker. Tagger and collaborators 
[ 31 ] have specifically argued that an outburst such as that shown in Figure [l] involves the destruction 
of a large-scale magnetic field via a reconnection event, during which the magnetic field eventually 
becomes weak enough that the magnetorotational instability can turn on. 

The simulations presented here begin with a value of a = 0.01, which is subsequently in- 
creased to a = 0.1 to mimic the effect of a jet ejection that removes the large-scale magnetic field 
in a particular region of the disk. In general, simulations of the magnetorotational instability sug- 



gest that a ~ 0.1 is an appropriate value when the instability is operating [15, m. The qualitative 



behavior discussed in the next section is not sensitive to the exact value of a chosen, however. 

It is important to point out the difference between this work and the wide variety of other 
models in the literature where a is allowed to vary in response to other disk parameters (e.g., 
[ p2[ , [TJ, |32|] for X-ray binary outbursts). The primary difference is that this work explores the 
effects of a rapid change in a due to an external event. As will be shown, the rapid nature of the 
change is a key ingredient that can drive an otherwise-stable disk into outburst. "Rapid" in this 
case means that a must change much faster than the disk is able to respond thermally, or vg S> ccc s , 
where vb is the speed at which the magnetic field is removed vertically from the disk (thereby 
leading to the change in a) and c s is the sound speed. This should not be difficult to achieve in a 
magnetically-dominated disk, where the Alfven speed is greater than the sound speed. 



3. Results 

Figure || shows the surface mass density evolution for a run in which the disk was initially held 
in a steady state with a = 0.01 for the first 100 seconds, and then a was increased to a value of 0.1 
in the inner part of the disk. The change in a drives the disk into outburst; the inner disk becomes 
nearly evacuated of matter, and a density wave propagates outward before eventually stalling when 
it reaches the region where a is still ~ 0.01. 

Figure ^ shows the evolution of the disk on a plot of surface mass density vs. temperature. 
The local energy balance curves for a = 0.01 and a = 0.1 are shown; the classic "S shape" of 
these curves is due to the interplay of radiation pressure and advective cooling We can see 
from this figure that although the disk is initially stable (and would remain so forever if a did not 
change), an outburst is triggered in the inner part of the disk simply because the change in a puts 
the disk in a region of parameter space (with respect to the new energy balance curve) where it can 
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Figure 2: The surface mass density evolution for a run in which the value of a is raised from 0.01 to 0.1 in 
the inner region of the disk (r < 400r g ; bottom panel). An outburst is initiated in the inner region, largely 
evacuating the disk and sending a wave of mass outwards, which eventually stalls when it reaches the region 
where a is still equal to 0.01. 
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Figure 3: The evolution on the local energy balance curve at two radii, when a is raised from 0.01 to 0.1 
in the inner region of the disk. The left panel shows the evolution at r = 22r g ; the energy balance curve for 
a = 0.01 is shown in black, and that for a = 0.1 is shown in blue. When a changes, the material at this 
location in the disk is suddenly in an unstable region of parameter space, and an outburst is triggered. At 
a radius of r = 327 'r g , meanwhile (right panel), the disk heats up when a changes, but it still remains on a 
stable part of the energy balance curve. Only when the density wave from the inner part of the disk reaches 
this radius (at ~ 140 seconds) is enough mass added to drive it into an outburst. 
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Figure 4: The outburst light curve for the run shown in Figures 2 and 3 (black), where the value of a is 
raised for r < 400r g , compared to the outburst light curve for a run in which the value of a is only raised for 
r < lOOr;, (red). The luminosity is plotted in units of the Eddington luminosity. 



heat up towards the advection-dominated upper branch. Figure [3] also shows that the outburst has a 
different character further out in the disk; there, the disk heats up when a is increased but does not 
immediately become unstable. Instead, the instability is triggered when the density wave arrives 
and the extra mass pushes this region of the disk "over the edge." 

The outburst light curve is shown in Figure ^ (black points). We see that the disk remains 
bright as long as the density wave is propagating through the disk and supplying mass to trigger new 
outbursts at large radii and simultaneously sustain the original outburst at the inner radii (cf. Figure 
When the density wave stalls, however, the inner radii stop receiving mass, and eventually the 
outburst dies down as the disk is able to cool. 

The density wave feature seen in these simulations is not new; it was observed in some of 
the earliest one-dimensional simulations of accretion disk outbursts (e.g., [17]) and in many others 
works since then. In those simulations, outbursts in an a = 0.1 disk were driven not by a change in 
a, but rather by an external accretion rate M that was high enough so that the equilibrium state of 
the disk would put it in the unstable portion of the local energy balance curve. In these simulations, 
the density wave stalls at a radius smaller than ~ 150 — 200r g , leading to an outburst of duration 



~ 20 - 30 seconds [|17|, gg, |30p. 

In the simulations I present here, there are several important differences. The timescale and 
radial extent of the outburst are not fixed, but rather depend on the size of the region in which a 
is changed. It is possible to obtain an outburst that is much longer than normal, as in the case 
discussed above, or much shorter, which will occur if we only allow a to change over a small 
region of the disk (as shown by the red points in Figure |]). Within a large range of radii over 
which a might be allowed to change, the density wave will always stall near the outer boundary 
of the region where a is altered (i.e., the outer boundary of the region where the "jet ejection" 
takes place). The reason that changing a can give longer outbursts than the original accretion disk 
simulations can be appreciated by looking at the right panel of Figure S by increasing a rapidly, 
we effectively obtain an a = 0.1 disk that has the surface mass density of an a = 0.01 disk. At 
any given radius, the disk is very close to being near the critical bend in the energy balance curve 
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at which it becomes unstable, and the density wave only needs to add a small amount of mass to 
trigger a new outburst at that radius. 

These results have important observational implications. The original simulations of Jl7| ] 
and others predicted continuously repeating outburst cycles, because the disks they studied were 
inherently unstable. GRS 1915+105, however (as well as other black hole candidates), frequently 
undergoes an outburst after a period of quiescence. The simulations presented here can reproduce 
this behavior easily, because it is the change in a that drives an otherwise stable disk into an 
outburst. Furthermore, if the value of M is such that the disk is stable for both the initial and final 
values of a, then we would expect to see a single outburst, with the disk returning to a stable 
state on the "new" energy balance curve after the outburst is complete. For a certain range of M, 
however, it is possible to have the disk be stable for the low value of a but unstable for the high 
value. In this case, the disk would be expected to undergo repeated oscillations following the initial 
outburst. Furthermore, these subsequent oscillations might be of the more "normal" variety found 
by [ jl"7| ] and others, because the disk at this point will be in a normal a = 0.1 state. Thus, one could 
imagine light curves with some resemblance to those shown in Figure [l| where the initial outbursts 
are of a different character and longer duration than the subsequent oscillations. Unfortunately, 
numerical instabilities associated with the inner disk boundary currently prevent investigation of 
the outbursts discussed here over more than one cycle, but I plan to pursue this in a future work. 

Figures |5| shows a run in which the value of a is only changed in the middle part of the disk, 
not the inner part. In this case, the change in a causes mass to build up at r ~ I00r g without 
penetrating to the extreme inner region. Examination of the local energy balance curves in this 
case shows that enough mass eventually builds up to trigger an outburst, but unlike in the case 
discussed above, this outburst is triggered by the existence of a gradient in a within the disk, not 
the speed with which a changes. 

The outburst light curve from this run is shown in Figure |6| The peak of the disk outburst 
is significantly delayed from the "jet ejection" (which occurred at ~ 100 seconds). This could be 
relevant for observations of GRS 1915+105, where infrared activity in the jet sometimes appears 
to precede significant X-ray activity from the disk (see Figure [T], or observations by Eikenberry and 
collaborators [||] for a more extreme example). These results can also be compared to the scenario 
investigated by [^0|], where a disturbance in the magnetic field at large radii led to an ejection as 
well as significant accretion into the inner disk — Figure ^ shows that when a is allowed to vary, 
the outburst does not necessarily have to reach the extreme inner disk in order to be observable. 

4. Conclusions 

A rapid increase in the strength of turbulence in an accretion disk (as might be associated with 
the removal of a large-scale magnetic field during a transient jet ejection) appears to be a promising 
way to explain how an otherwise-steady disk can be driven into outburst, as well as to explain 
some of the details of the outburst light curves seen in GRS 1915+105 and other X-ray binaries. 
Future work should involve a more rigorous comparison of the model to observations, as well as 
comparison to X-ray spectra. Overall, the period of activity in an accretion disk immediately after 
a jet has been ejected may be a promising place to test theoretical models, because the physics of 
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Figure 5: The surface mass density evolution for a run in which the value of a is raised from 0.01 to 0.1 
in the middle part of the disk (cf. Figure 2). In this case, mass builds up around the boundary of the region 
where a is changed and penetrates it slightly, but never reaches the inner part of the disk. Also, unlike Figure 
2, there is no prominent density wave initiated when a changes. 
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Figure 6: The outburst light curve for the run shown in Figure 5. The outburst peak is significantly delayed 
from the time (t = 100 seconds) at which a was changed. 
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the disk in this phase is well-understood but the behavior of the disk may still be influenced by 
fingerprints of the jet. 
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